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We study the single transverse spin asymmetries in the single inclusive particle production within 
the framework of the generalized parton model (GPM). By carefully analyzing the initial- and 
final-state interactions, we include the process-dependence of the Sivers functions into the GPM 
_ ^ formalism. The modified GPM formalism has a close connection with the coUinear twist-3 approach. 

. Within the new formalism, we make predictions for inclusive tv^ and direct photon productions at 

RHIC energies. We find the predictions are opposite to those in the conventional GPM approach. 
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I. INTRODUCTION 

rS i' Single transverse-spin asymmetries (SSAs) in both high energy lepton-hadron and hadronic scattering processes have 
JL , attracted considerable attention from both experimental and theoretical communities over the years Generally, 
defined as = icr{S±) — a{—S±))/{a{S±) + a{—S±_)), the ratio of the diff'erence and the sum of the cross sections 
' when the hadron's spin vector Sj_ is flipped, SSAs have been consistently observed in various experiments at different 
'— ' . collision energies 0-01 • 

Much theoretical progress has been achieved in the recent years. An important realization is the crucial role of the 
initial- and final-state interactions between the struck parton and the target remnant Q , which provide the necessary 
phases that leads to the non-vanishing SSAs. These interactions can be accounted for by including the appropriate 
, color gauge links in the gauge invariant transverse momentum dependent (TMD) parton distribution functions (PDFs) 
0^ ' 043 ■ An important example is the quark Sivers function 9] , which represents the distribution of unpolarized quarks 
in a transversely polarized nucleon, through a correlation between the quark's transverse momentum and the nucleon 
• polarization vector. They are believed to be (partially) responsible for the SSAs observed in the experiments. 
' The details of the initial- and final-state interactions depend on the scattering process, thus the form of the gauge 
. link in the Sivers function is process dependent [lo| . As a result, the Sivers function itself is non- universal. For 
example, it is the difference between the final-state interactions (FSIs) in semi-inclusive deep inelastic scattering 
(SIDIS) and the initial-state interactions (ISIs) in Drell-Yan (DY) process in pp collision that leads to an opposite 
sign in the Sivers function probed in these two processes j6l. Isl. [lH . For hadron production in pp collisions, typically 
the Sivers function has a more complicated relation relative to those probed in SIDIS and DY processes [lOj; that 
is, there arc only FSIs (ISIs) in the SIDIS (DY) process, while both ISIs and FSIs exist for single inclusive particle 
" " ' production. 

The SSAs for inclusive single particle production in hadronic collisions are among the earliest processes studied 
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in experiments, starting from the fixed-target experiments in 1980s [12|. Recently the experiments at Relativistic 
Heavy Ion Collider (RHIC) have also measured the SSAs of inclusive hadron production in pp collisions over a wide 
range of energies (^]. Theoretically a QCD collinear factorization formalism at next-to- leading-power (twist-3) has 
been developed and been used in the phenomenological studies p"3l - fl6j . Alternatively, a more phenomenological 
approach has also been formulated in the context of generalized parton model (GPM) (T7l-[T9|. with the inclusion 
of spin and transverse momentum effects. In this approach TMD factorization is assumed as a reasonable starting 
point [13]; at the same time, the leading twist TMD distributions (Sivers functions) are assumed to be universal 
(process-independent), thus the same as those in SIDIS process [2ll[2^. 

In this paper we formulate the SSAs in inclusive single particle production within the framework of the GPM 
approach. However, instead of using a process-independent Sivers function, we will carefully examine the initial- and 
final-state interaction effects, and determine the process-dependent Sivers function. Further we find one can shift the 
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process-dependence of the Sivers function to the squared hard partonic scattering amphtude under one-gluon exchange 
approximation, and these modified hard parts are very similar in form as those in the twist-3 colhnear approach [l5[ 
in terms of Mandelstam variables s,t,u (as we will demonstrate). This suggests a close connection between this 
modified GPM formalism and the twist-3 approach. However, it is important to mention that Mandelstam variables 
s,t,u are themselves a function of partonic intrinsic transverse momentum in the GPM approach. We comment 
on these issues at the end of Section [TTl where we also show the modified GPM formalism can reproduce the twist-3 
collinear factorization formalism in the leading order expansion in intrinsic transverse momentum kx (for contributions 
coming from initial and final state interactions, where the latter is equivalent up to a prefactor). The rest of the paper 
is organized as follows: In Sec. |lTl we introduce the GPM approach, demonstrate how to formulate the ISI and 
FSI effects, and discuss the connection to the twist-3 collinear factorization approach. In Sec. Illli we estimate the 
asymmetry for inclusive pion and direct photon production at RHIC energy, and compare our predictions with those 
from the conventional GPM approach. We conclude our paper in Sec. IIVI 



II. INITIAL- AND FINAL-STATE INTERACTIONS IN SINGLE INCLUSIVE PARTICLE 

PRODUCTION 

In this section, we introduce the basic ideas and assumptions of the GPM approach. Then we discuss how to 
formulate the initial- and final-state interactions for single inclusive particle production. Within the same framework 
of GPM approach, we thus derive a new formalism for the SSAs of single inclusive particle production, with the 
process-dependence of the Sivers function taken into account. 



A. Generalized Parton Model 



The generalized parton model was introduced by Feynman and collaborators i23| as a generalization of the usual 
collinear pQCD approach. It was adapted and used to describe the SSAs for inclusive particle production |T7l - [l9| . 
which has had considerable phenomenological success [isj . According to this approach, for the inclusive production 
of large P^t hadrons (or photons), A'^{Pa) + B{Pb) — >■ h{Ph) + X, the differential cross section can be written as 
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where S = {Pa + Pb)^, fa/A^i^a^kaT) is the TMD parton distribution functions with kaT the intrinsic transverse 
momentum of parton a with respect to the light-cone direction of hadron A, and Di^/^{zc) is the fragmentation 
function. Since we will only consider the SSAs generated from the parton distribution functions in this paper, we 
have neglected the /cT-dependence in the fragmentation function. -ff^_j.j,(s, t, u) is the hard part coefficients with s, t, u 
the usual partonic Mandelstam variables. Eq. ([1]) can also be used to describe direct photon production, in which one 
replaces the fragmentation function Di^/^{zc) by d{zc — 1), and by aemO-s- 

To clearly specify the kinematics, we consider the center-of-mass frame of the two initial hadrons, in which one 
has P% = ^fSj2n^ and = .Jsj2n>' , with fi^ = [l+,0-,0±] and 71^ = [0+,l",0±] in light-cone components. 
For future convenience we also define the hadronic Mandelstam invariants, T = [Pa — Ph)^ and U = {Pb — PhY ■ 
Additionally, the momenta of the partons in the partonic process a{pa) + b{pb) — >■ c{pc) + d{pd) can be written as 
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where the momentum of parton c is related to the final hadron as: pc = Ph/zc- 

To study the SSAs, the PDFs fa/A^{^a, kar) in the transversely polarized hadron A can be expanded as [T7l-[2C 

fa/A^{Xa,kaT) = f a / a{x a , klj,) + f a , klj,) - 
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where Sa is the transverse polarization vector, M is the mass of hadron A, fa/A{xa, k'^rp) is the spin-averaged PDFs, 
and f^{xa, fc^r) Sivers functions. Thus in GPM approach, the spin-averaged differential cross section is given 

by Eq. ([T]) with fa/A^ i^a, kar) replaced by fa/A{^a, fcaT); while the spin-dependent cross section is given by 
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and the SSA is given by the ratio, 



„ dAa I „ da 



As stated in the introduction, there are two assumptions in the GPM approach: one is that the spin-averaged 
and spin-dependent differential cross sections can be factorized in terms of TMD PDFs as in Eqs. ([T]) and 
and the other one is that the Sivers functions is assumed to be universal and equal to those in SIDIS process, 
f^{xa, k^x) ~ /it* '^^^^^(-^a, fcar)- ^^^^ paper we continue to work within the framework of the GPM approach, 
in other words, we will assume the TMD factorization is a reasonable phenomenological starting point. However, at 
the same time, we will take into account the initial- and final-state interactions. Since both ISIs and FSIs contribute 
for single inclusive particle production, in principle the Sivers functions in inclusive particle production in hadronic 
collisions should be different from those probed in SIDIS process. We thus need to carefully analyze these ISIs and 
FSIs for all the partonic scattering processes relevant to single inclusive particle production to determine the proper 
Sivers functions to be used in the formalism. In other words, this new formalism will be 
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in which a process-dependent Sivers function denoted as f^''^''^'^'^{xa,k^rp) is used rather than that from SIDIS 
y.i^,siDis^^^^ ^2^-j conventional GPM approach. 



B. Initial- and final-state interactions 



In this subsection, we will discuss how to formulate the initial- and final-state interactions. The crucial point is that 
the existence of the Sivers function in the polarized nucleon relies on the initial- and final-state interactions between 
the struck parton and the spectators from the polarized nucleon through the gluon exchange. Thus by analyzing 
these interactions, one can determine the process dependent Sivers function /-(^'"''"'^'^''(xa, fc^T) to be used for the 
corresponding partonic scattering ab — >■ cd. We start with the classic examples: the final-state interaction in SIDIS, 
and the initial-state interaction for DY process. To the leading order (one-gluon exchange), they are shown in Fig. [1] 
For the SIDIS process e{£) +p{Pa, St) — >■ e{£') + h + X with = —q^ — — (£' — under the eikonal approximation. 




FIG. 1: Final-state interaction in SIDIS (left) and initial-state interaction in DY (right) processes. 



the final-state interaction (as in Fig. [TJleft)) leads to 



-i{pc){~i9)l-T- 
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u{pc) 
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(7) 



where the gamma matrix 7" appears because of the interaction with a longitudinal polarized gluon (~ ^^), and a 
is the color index for this gluon. The eikonal part (the term in the bracket) is the first order contribution of the 
gauge link (in an expansion of the coupling g) in the definition of a gauge-invariant TMD PDFs in SIDIS process, see 
Fig. [2][a). The imaginary part of the eikonal propagator l/(— fc+ -|- ie) provides the necessary phase for the SSAs. 
On the other hand, for DY process, the initial-state interaction (as in Fig. [Ijright)) leads to 
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FIG. 2: Sivers function in SIDIS process in the first non-trivial order (one-gluon exchange). 



which has the same real part and opposite imaginary part compared to SIDIS process. This leads to the fact that the 
spin-averaged TMD PDFs are the same, while the Sivers function will be opposite in SIDIS and DY processes. This 
conclusion can be generalized to all order, and has been proven to be true using parity and time-reversal invariant 
arguments [1, Q. 
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FIG. 3: Initial- and final-state interactions in qq' qq': (a) initial-state interaction, (b) final-state interaction, (c) and (d) the 
final-state interactions for the unobserved particle. 



Now let us turn to the case for inclusive single particle production in hadronic collisions, in which 2 — )> 2 partonic 
scattering is the leading order contribution, where both initial- and final-state interactions contribute. We will 
start with a simple example: qq' — > qq' . Here the initial-quark q is from the polarized nucleon, and the final-quark q 
fragments to the final-state hadron. The one-gluon exchange approximation for the initial- and final-state interactions 
are shown in Fig. [3l Under the eikonal approximation, for ISI Fig.[3][a), 
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Likewise, for the FSI Fig. [S^b), we have 
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Thus both interactions contribute to the phase — i7r(5(fc"'"), which is the same as in the SIDIS process as in Eq. ([T]). 
However, they will have different color flow. To extract the extra color factors for Fig. BTa) and (b) as compared to 
the usual qq' -> qq' without gluon attachments, we resort to the method developed in jl4 [isl [26j . We obtain the 
color factors C/ [Cp^) for initial (final)-state interaction 
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while the color factors for unpolarized cross section is given by 
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In other words, the Sivers function in qq' qq' should be the one as shown in Fig. |4l which comes from the sum of 
the ISIs and FSIs with the corresponding color factors C/ and Cp^ respectively. Thus by comparing the imaginary 
part of the eikonal propagators in Eq. ([7]) for SIDIS and those in Eqs. © and for ISI and FSI for qq' qq', we 
immediately find the Sivers function probed in qq' — !■ qq' process is related to those in SIDIS as follows 



p±a,qq'^qq' 
JIT 



Ci 



y._La. SIDIS 
JIT 



(13) 



5 



P. 




(a) (b) 
FIG. 4: Sivers function in qq' — >■ qq' from ISIs and FSIs, with tiie corresponding color factors Ci and Cf^ respectively. 



Thus in the GPM model, using the process dependent Sivers function, one should replace 
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by the following form 
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where hqqi-^qqi is the partonic cross section without color factors included. For qq' — >■ qq' , one has 
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Alternatively one can use /j^'^^^^^ for the single inclusive particle production while accounting for the process- 
dependence of the Sivers function, by shifting the process-dependence to the hard parts. In other words, instead of 

qq'^qq' 



using H^Li^qq, in Eq. ^ for the spin-dependent cross section, one should use 



where 
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are the corresponding hard parts related to initial- and final-state interactions, respectively. 

There are many other partonic processes contributing to the single inclusive particle production. Similar to the 
analysis in qq' — qq' ^ one needs to analyze each individual Feynman diagram accordingly, carefully moving the extra 



factors (process-dependence) from the corresponding Sivers function to the hard parts, thus obtaining H^^^^ and 
^ab^cd ^'^^ every channel. The modified formalism will be given in the next subsection. 

There are some comments to our results presented to this point: in particular those displayed in Fig. 31 It looks 
like Figs. ISja), (b) can be factorized into a convolution of Sivers function and a hard part function as shown in 
Fig. m However, this is not a TMD factorization in the strict sense. Currently TMD factorization theorems have 
been established for both SIDIS and DY processes [H, [2^. To the order we are studying, this means, the one-gluon 
exchange diagram for SIDIS in Fig. [T]can be factorized into a convolution of a Sivers function f^'^^^^^{xa, ^ot) ^'^d 
a hard part function H{Q), as shown in Fig. [2j Here all the soft physics (those depending on fear) has been absorbed 
into the Sivers function f^'^^^^^{xa,k'^x)' ^^'^ hard part function H{Q) only depends on the hard scale Q, not 
kaT- On the other hand, for qq' — >■ qq' ^ we write the corresponding diagram Fig. [2Ia) into a similar form: a product 

of a Sivers function /j^'*' (xa^ k'^r) ^^'^ ^ hard part function Hqqi ^qqi[s, i, u), as shown in Fig.|4l But as we will 
comment later, besides the kaT dependence from the Sivers function, one will also need to keep the kaT dependence 
in the hard part functions Hqqi^qqi, without which the SSAs will vanish in both the conventional GPM and this 
modified GPM formalism. Even though this is not a TMD factorization, one hopes this formalism is a reasonable 
approximation. There are two reasons to suggest this might be the case. First of all, from phenomenological point of 
view, this formalism had some success [ll|. Secondly, as we will show in Section |lI]D this formalism has a connection 
with the well-established collinear twist-3 approach [l^. In this respect, our identification of the color factors with 
the hard cross-sections is reminiscent of the results of the twist 3 approach (see in particular [ill). Indeed we will 
see that upon calculating all partonic processes that contribute from each channel, they have the same form in terms 
of Mandelstam variables s, i, ii, as compared to those in the twist-3 collinear factorization approach [TBj (up to a 
prefactor associated with final state interactions). 
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To close this subsection, we want to point out the fohowing important fact: the interaction with the unobserved 
particle (the quark for qq' qq') vanishes after summing different cut diagrams 0, To see this clearly, 

we have for Figs. Wic) and [3l^d) 

l^^5{pl)^^in5{{pa-kf)5{pl), and -±-5{{p^ - kf) ^ +iiT5({pd - kf)5{pl). (19) 

[Pd-ky+ie Pd^'^^ 

respectively. Since the remaining parts of the scattering amplitudes for these two diagrams are exactly the same 
except for the above pole contributions which are opposite to each other, the contribution from the unobserved 
particle vanishes. This could also be used to explain why the inclusive DIS process, the SSA vanishes. As shown in 
Fig. [T] (left), we don't observe the final-state quark for the inclusive DIS process, thus the contribution from the cut 
to the left and to the right will cancel which results in a vanishing asymmetry. 

We want to emphasize that the above analysis holds true only under one-gluon exchange approximation. Going 
beyond one-gluon exchange, the Sivers functions are typically more complicated, there seems no simple relation (as 
extra color factors) to those in the SIDIS process psj . 



C. Single inclusive hadron production 



Now after carefully taking into account both initial- and final-state interactions, the more appropriate GPM for- 
malism for spin-dependent cross section should be written as 
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where we have a new hard part function H^^'l^^ instead of H^j^^^ used in the conventional GPM approach. Here the 
process dependence in the Sivers function has been absorbed into -ff^"^^' which can be written as 
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where H^b%^ and H^'^^ are associated with initial- and final-state interactions, respectively. The contributions for 
the various contributing partonic subprocesses are given by 
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We also calculate the corresponding hard part functions for direct photon production, and they are given by 



H, 



Inc 

qg^jq 



-H. 



Inc 

qg^iq 



rrlnc _ TT^^^^ 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



Here again we note that all these hard part functions have the same form in terms of Mandelstam variables s, i, u, 
compared to those in the twist-3 collinear factorization approach [l^: H^"^^ and H^^^~J have the same functional 



form as the corresponding ones H^]^~ ~ and H^^~ ~ (defined below) in the twist-3 collinear factorization for- 
malism, respectively. However, there are two differences in the formalisms. First, in the twist-3 collinear approach, 
the hard part functions are given by 
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twist 3 F associated with final state 



i.e., there is an extra factor (1 -I- u/i) accompanying the hard part functions i/^'^ 
interactions. However, in our modified GPM formalism as in Eq. ([2T|) . there is no such factor. This difference can be 
traced back to the eikonal approximation we are using, see, e.g., Eq. (jlOp . where we only keep the pole contribution 
— fc+ + ie in the denominator under this approximation. However, there is an extra term linear in k± (cx Pc ■ k±) 
which exists in the twist-3 collinear factorization formalism. This leads to the extra factor (1 -|- u/t) for the final-state 
interaction contribution (for details, see Ref. [Tsj). Second, in the twist-3 collinear factorization approach, all the 
parton momenta are collinear to the corresponding hadrons, thus s, t, u does not depend on the parton intrinsic 
transverse momentum. On the other hand, in the GPM approach the parton momenta involve intrinsic transverse 
momentum, thus s, t, u all depend on the the parton transverse momentum, kaT and ki,T. In fact, because of the 
existence of the linear fcoT-dependence in efc-TSAnft^ one has to keep another linear /caT-dependence from the rest of 
the integrand in Eq. (|20p . otherwise the integral over cPkaT vanishes. In other words, it is the linear in kaT term in the 
hard part functions H^'l'i^^{s, t, u) and d{s + t + t) that contributes to the asymmetry. Even with these two differences, 
the similarities in terms of s, t, u suggest that there are close connections between our modified GPM formalism and 
the twist-3 collinear factorization approach. We explore this potential connection in the next subsection. 



D. Connection to the twist-3 collinear factorization formalism 



As pointed out in the last subsection, it is the linear in kaT dependence from the rest of the integral in Eq. ((20|) 
that contributes to the asymmetry. We thus make an expansion and keep only the linear in kaT terms. Wc will show 
that the leading term in this expansion has a close connection to the twist-3 collinear factorization formalism. 
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We start by specifying the partonic kinematics. Keeping the hnear in kaT terms and dropping ah the feb^- dependence 
we have ~ XaP^ + kaT and p'^ w xi,Pg, thus 
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and X = —XbU / {zcXbS + T) is independent of kaT- Now performing the integrate over Xa in Eq. (j20p and using the 
i5-function we get, 



Eb 



dAa 



S 



/ d kaT J7 — /it {Xa,k^j 



a.b,c 

dz, 



M Xa 



dxi 

Xb 



fb/Bixb) 



XbS + T/zc 



(38) 



After replacing above, one has 



s = s - -2F/iT • kaT/ Zc, i = i+ -2PhT ■ kaT/ Zc 



(39) 



where s — xxbS, t = xT / z^-, u — XbU / z^ and they are all independent of kaT- Note s + i + u = implies s + i + u = 0. 
Now besides the g'^^TSAnn^ ^j^g linear in kaT contributions in Eq. (p8| can come from, either (a) aia-dependence in 
jj_a,siDiSj-_^^^ ^2_^-j^ or (b) the I- and i-dependence in i7^"l^^(s, i, m). This is because Xq, s, and i are the only terms 
in Eq. psp which depend linearly in kaT- We now make kaT expansion one by one. First for contribution (a), since 



dXa _ 

dk'^rp ZcXbS + T' 



(40) 



to the linear term in kaT, we have 



Eh 



d^Ph 



- f E / 



a,b,c 
2 



(zkaTSAnn 



M 

rinc 



2P, 



hTa 



aT 



ZcXbS + T dXa 
1 



.±a,SIDIS/ ,2 
J IT \Xa,i^al 



dXb , , . 

Jb/B[Xb) 

Xb 



(41) 



where we have dropped all kaT dependence in H^'z^^, thus replacing the fcoT-dependent s, t, ii by the fcaT-independent 
s, i, u in H^^^_^^. Then using 



/j2; 1 13 J a /.±a,SIDIS/ 72 \ ^ f j2 1 0a iT" 1 2 j- J-a,SIDIS ^ ;2 \ 

d kaTk'aTkaTflT {Xa,k^j.) ^ -- d kaT d'^ I ^aT I fiT {Xa,kaT), 

and the relation between the Sivers function and the Efremov-Teryaev-Qiu-Sterman function Ta,Fix,x) 



Ta,Fix,x) 



1 

M 



d'kaT\kaT\'f^T'''''''{x,k'iT), 



(42) 



(43) 



one can rewrite Eq. (|41l) as 

a,b,c 



ZcU X 



Ta,F{x, X) - X — Ta,F{x, X) 

dx 



^fb/Bixb)H':r-.Mi,u)- 

Xb 



XbS + T/zc 
(44) 

We observe that this form is the same as that in the twist-3 collinear factorization approach. In particular, note that 
there is no fcaT-dependence in the hard part functions i?^^!^^. The difference to the twist-3 collinear factorization 
formalism [l5| (as mentioned above) is the extra factor (1 + u/i) accompanying the hard part functions associated 
with final-state interactions, see Eqs. (I^Tj) and ([55)1 . 
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However, in our modified GPM formalism, we liave anotlier contribution from (b), due to the /caT-dependence from 
H^%^{s,t,u) in Eq. (|38| . Let's now study this contribution (b). As is explicit in Eq. (p9|) u is independent of kaT 
while both s and i depend on kaT- Since s + t + u = Q, one could then set t — ~s — u in -ff^b^-c then expand only 
s in kaT- That is, 



d 



dk 



rrinc 



(s, i, u) 



aT 



ds d 
dk^di' 



H. 



Inc 



U, U ' 



2sPhTa d 



ds 



(45) 



Then we have the contribution (b) 



E, 



d^Ph ~ S 



J2 J ^Dh/c{Zc) 

a,b,c 



r Ta^F[X,X) / h/B[Xb) 

ZcU X J Xb 



d 

^S-^Hll%^{S,-s ~ u,u) 



1 



XbS + T/zc' 
(46) 

Thus to the leading order (linear in kaT terms) , the spin-dependent cross section in our modified GPM formalism can 
be written as 



^ dAa ^ dAcr(°) 



d^Pb 



d^Ph 



Eh 



'dm^ 



(47) 



with the contributions (a) and (b) given by Eqs. (|44p and respectively. The term (a) almost reproduces the twist-3 
collinear factorization formalism in Ref. [15] modular the extra factor (l + u/t) associated with final state interactions, 
for which the origin of the difference is understood in last subsection. On the other hand, for the extra term (b), 
theoretically how to interpret this "mismatch" and why the term (b) does not appear in the usual twist-3 collinear 
factorization formalism deserves further investigation Here it is important to note, from the phenomenological 
perspective, as already shown in fisj . the derivative of the correlation function Ta^rix, x) is the dominant contribution 
to the SSAs, thus we expect the term (b), which contains no derivative, to play a less important role in generating 
the SSAs compared with term (a). In other words, even though this modified GPM has an extra piece compared with 
the well-known twist-3 collinear factorization formalism, phenomenologically (numerically) this formalism could give 
a good approximation to the SSAs. This remains to be confirmed [29] because there is still a difference in term (a) 
on the extra factor (1 -I- ii/t) associated with the final state interactions between the twist-3 collinear factorization 
approach and our modified GPM formalism. If this were the case, it will provide further support to the modified 
GPM approach to the SSAs. 

To close this section, we want to emphasize that the contribution calculated in Ref. jT^ only comes from the 
so-called soft-gluon-pole (SGP) in the twist-3 collinear factorization approach. However, there are also contributions 
from so-called soft-fermon-pole (SFP) js^j- Even though our modified GPM formalism might capture the main feature 
of SGP contributions, it seems unlikely to reproduce the SFP contributions. In this respect the twist-3 formalism is 
"internally complete" in the sense that the collinear factorization is expected to hold for this formalism [3l| . Finally, 
while TMD factorization is assumed in both GPM and our modified GPM formalisms, it is likely not to hold in 
these processes J2§\. However, the extent to which it is broken is not known numerically. Thus, calculations within 
(modified) GPM formalisms should bear this in mind and thus be used with extra care. 



III. NUMERICAL ESTIMATE OF THE SSAS 



In this section, we will estimate the SSAs for single inclusive hadron and direct photon production in pp collisions 
at RHIC energy by using our modified GPM formalism in Eq. I^U^ . We will compare our results with those calculated 
from the conventional GPM formalism as in Eq. 

To calculate the spin-avera ged cross section, we use GRV98 LO parton distribution functions [32?] along with a 
Gaussian- type fcr-dependence |2TI. [22| . The hard part functions for different partonic scattering channels are available 
in the literature [H, [H, 113] • For the spin-dependent cross section, we use the latest Sivers functions from which 
are extracted from the recent SIDIS experiments. To consistently use this set of Sivers function, we will use DSS 
fragmentation function [35| . For the numerical predictions below, we work in a frame in which the polarized hadron 
moves in the -|-z-direction, choosing S±, Phi_ along y- and a;-directions, respectively, where all the relevant distribution 
functions and fragmentation functions evaluated at the scale Ph± [iTj . 

In Fig.O we plot the Ajv as a function oi xp for inclusive 7r° (left) and direct photon (right) production at rapidity 
y = 3.3 for RHIC energy y/s = 200 GeV. The estimates using the conventional GPM formalism in Eq. (j4|) are shown 
as dashed lines, while those using our modified GPM formalism in Eq. ([20)1 are shown as solid lines. One immediately 
see that for both inclusive tt" and direct photon, A^r change signs compare to the conventional GPM formalism. For 
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FIG. 5: An for inclusive particle production as a function of a:F at RHIC energy ^ = 200 GeV: p^p — ^ 7r° + X (left) and 
p^p — >■ 7 + X (right). The dashed curves are for the conventional GPM calculation, and the solid curves are for our modified 
GPM calculation. We have used the latest Sivers function from 22?|, and DSS fragmentation function [35l ]. 



TT^ , the conventional GPM predicts a negative asymmetry (though very small from this set of Sivers functions) , while 
the modified GPM formalism predicts a positive asymmetry. On the other hand, for direct photon, conventional GPM 
formalism predicts a positive asymmetry, while modified GPM formalism predicts that the asymmetry is negative, 
which is consistent with the predictions from twist-3 collinear factorization approach fl^. This can also be easily 
understood as follows. In the conventional GPM approach, one use in the calculation of the spin-dependent cross 
section. For direct photon production, the dominant channel comes from qg — > jq, with jl5l , |33| 



ttu _ _;_ 



t s 

s i 



while the hard part in the modified GPM formalism is given by 



H, 



Inc 



Nr. 



t s 

s i 



(48) 



(49) 



This introduces an extra color factor —N^/{N^ — 1), thus opposite to the conventional GPM formalism. This prediction 
comes from the process-dependence of the Sivers functions, and has the same origin as in the photon-|-jet calculation 
[36j . On the other hand, for the inclusive production, the dominant channel comes from qg qg, particularly in 
the forward direction, one has 



rrlnc rrlnc — I i rrinc — F 

qg^qg ^ qg^qg qg^qg 



2s' 



1 



2s' 



2 s' 



2{N^ - 1) £2 N2 _ I £2 



m -It'' 



where we have used that in the forward direction, i is small, while u ■ 



^ _ N',~l 
qg-^qg 2N' 



s 
u 



t' 



whereas Il5l |3 
2s' 



(50) 



(51) 



We thus also see the sign is reversed in our modified GPM formalism compared with the conventional GPM approach. 

We observe that the XF-dependence in both modified and conventional GPM formalisms are different from those 
observed in the RHIC experiments where larger asymmetries have been observed in the forward direction (large 
xf) H- Of course, in order to have a comparison with the experimental data for inclusive hadron production at RHIG 
experiments, one must include both Sivers (as studied in this paper) and Collins effects js^]. The latter describes a 
transversely polarized quark jet fragmenting into an unpolarized hadron, whose transverse momentum relative to the 
jet axis correlates with the transverse polarization vector of the fragmenting quark. This latter correlation can also 
generate the transverse spin asymmetry (which is not studied here). Currently attempts at global fitting with both 
SIDIS and pp experimental data are ongoing [l^. We encourage the use of the modified GPM formalism in such 
a global analysis, to study the effect of the associated ISIs and FSIs (process-dependence of the Sivers functions). 
We also emphasize '36'] that there is only Sivers contribution in direct photon production. Since the modified and 
conventional GPM predict opposite asymmetries, direct photon production presents a favorable opportunity to test 
the process dependence of the Sivers function, or the effect of the associated ISIs. 
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IV. SUMMARY 



In this paper, we have studied the single transverse spin asymmetries in the single inclusive particle production 
in hadronic collisions. We point out the Sivers functions in such processes are generally different from those probed 
in the SIDIS process because of different initial- and final-state interactions. By carefully taking into account the 
process-dependence in the Sivers functions (under one-gluon exchange approximation), we derive a new formalism 
within the framework of GPM approach. We find this formalism has close connections with the coUinear twist-3 
approach. With our modified GPM formalism, we make predictions for the inclusive tt'^ and direct photon production 
in pp collisions at RHIC energies. We find that the asymmetries predicted from the modified GPM formalism are 
opposite to those in the conventional GPM approach. This sign difference comes from the color gauge interaction, 
which has the same origin as the sign change for Sivers functions between SIDIS and DY processes. Our predictions 
about the sign are consistent with those from the twist-3 collinear factorization approach. We encourage a global 
analysis of both SIDIS and pp experimental data using this modified GPM formalism. 
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